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Notation 


DK, Fluid film radial viscous damping and stiffness, 
respectively 

i}, A Rotor modal external damping, stiffness, and 

mass, respectively 

Fluid dynamic force in coordinates rotating at 

the rave AR. 

yT 


F. =F, +jF, 


i 
f Time 
209 =x +A Rotor lateral displacement with its horizontal 
Ce) and wertical (7) components in Stationary 
coordinates. 
af = Ey +in(f) Rotor lateral displacement with its orthogonal 
components £, in the coordinates rotating at 
the rate Ag 
Fluid circumferential average velocity ratio 
Rotor rotative speed 
a 
dt 
a” 


ar 





models of rotors, this new rotor/fluid system model opens 
new paths for the prediction of rotor stable rotation and the 
prevention of damaging self-excited lateral fluid whirl and 
whip vibrations [Ref. 4, 5]. 


Fluid force model 


The physical situation of the rotor/bearing or rotor/seal 
system can be modeled by two coaxial cylinders, one softly 
supported and rotating inside the other, stationary cylinder. 
The relatively small clearance between these cylinders is filled 
with fluid (gases included). Due to the rotation of the inner 
cylinder’, and, due to friction, the fluid within the clearance 
becomes involved in circumferential motion. 


Fluid dynamic forces in the clearance can be modeled in the 
form of a dashpot and a spring, as is illustrated in Figure 1 (for 
the moment, we will ignore the fluid inertia). The fluid film- 
related dashpot and spring are not, however, stationary, but 
rotate at the rate AQ, a fraction of the rotative speed & (capital 
omega). The fraction parameter, A (lambda), has a physical 
meaning; it is the fluid circumferential average velocity ratio. 
Figure 2 illustrates the fluid velocity profile and the average 
circumferential angular velocity. 

Introduce two systems of reference coordinates: X and Y are 
stationary coordinates, and the coordinates & (xi) and 7 (eta) 
are coordinates that rotate together with the dashpot and 
spring system, at the fluid average circumferential angular 
velocity AQ. The fluid force components can be simply 
expressed in these rotating coordinates as: 

k 

F =DE+K E (1) 

F,=Dq+K,n 
or 


Fe=FetjF,=DE+K,¢, f=V-T, UO =E +n (2) 


where? is time; F, is the fluid force in the rotating coordinates 
En; D is fluid viscous damping; K, is fluid film radial 


his article summarizes the results of experimental 
modal testing of rotors in fluid-lubricated bearings and 
seals which was performed at Bently Rotor Dynamics 
Research Corporation (BRDRC) during the last ten years (Ref. 
1]. A new and efficient model of the fluid forces in seals and 
lightly-loaded, fully lubricated bearings emerged from these 
results [Ref. 2, 3]. Used together with multi-mode, modal 





Figure 1 
Fluid force generated inside the clearance: the spring and 
dashpot rotate at the rotative speed Af. 


1 Avery similar situation occurs when the outer cylinder mags, and the inner one is stationary, or when both cylinders rotate at different speeds, such as in the case of partially engaged 
hydraulic clutches, 
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stiffness; G(zeta) is the rotor lateral displacement using 


complex rotating coordinates and « is pr 
f 

Traditionally, the dashpot damping force is proportional to 
the velocity of the dashpot end at the rotor surface, and the 
stiffness force is proportional to the spring end displacement. 
Note, however, that the spring and dashpot ends are not 
attached, either to the bearing or to the rotor surfaces, but 
slide with the velocity ANR +c) on the bearing and (1-A) QR 
on the rotor, where R is the rotor radius and c is the radial 
clearance, 


Rotor model 

Consider the rotor itself (Figure 3). The isotropic, balanced 
rotor is represented by its first lateral mode modal parameters: 
stiffness K , mass M , and external damping, D, Omitting the 
fluid force for a moment, the rotor model — its equation of 
motion, a balance of forces expressed in the stationary coordi- 
nates XY, is: 

Mx+D iX+Kx=0 

My+D.jp+Ky=0 
or 

Mz+D,z+Kz=0 z=x+jy (3) 
where z(t) =x(t) +t) is the rotor lateral complex deflec- 
tionand = is —— 


ar 

Eq. (3) can be transformed into rotating coordinates En, 
using the following transformation algorithm: 

x= EcosAgtnsinAm 

y= Esini + cosà 
or, further applying Euler's identity e’*=cose +/sin a: 

z=fe"™; zaxthy, 2=E+/n (4) 
: For Eq. (3) transformation, the velocity? and acceleration 
7 are required: 
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Figure 2 
Fluid velocity profile inside the bearing or seal clearance. 
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. dz Ls 
Z=—=(F +janter™ 
at 
bè dz +? . 
z= = (2+ 2jAnz-A*n* Ze" ™ 
dt” 


Eq. (3) in the coordinates rotating at the rate AQ will have, 
therefore, the following form: 


MÈ + 2jAME-A707Z) + D AÈ +JAQZ) +KZ=0 (5) 
The first term in Eq. (5) is composed of three parts; the 


relative inertia force, the Coriolis inertia force, and the cen- 
tripetal inertia force. 


Model of a rotor with fluid interaction 

Now Eq. (5) can be supplemented by the fluid force (2), as 
all forces are expressed now in the same coordinate system. In 
rotating coordinates £n, the rotor/fluid model will have the 
following form: 

M(Z + 2fAQZ-A707Z) + D(Z+jANTZ) +KT+DZ4+K,2 =0 

(6) 

To transform the rotor/fluid model back to the stationary 
coordinates, the inverse transformation (4) must be used: 

T (7) 

The rotor/fluid model in the stationary coordinates has its 
final form: 

MZ+D,Ž+Kz+D(JAQZ) +K,z=0 (8) 

The last three terms of Eq. (8) represent the rotating fuid 
force. 

When the lateral complex coordinate z is split back to the 
horizontal x and vertical y displacements, Eq. (8) will have the 
following familiar form: 

Mx +D X +Kx+ D(x + AQ) +K,x=0 

MV +D y+ Ky + D(p-aQx) +K, 7 =0 
or, using the matrix notation: 


(9) 


K+K,  aad|lx 
+ + =0) 
l-20 K+K, Ily 


w 
+ 
I 
= 
AT] 


2 
S 
+ 
=, 

iz 





(10) 





Rotor model (bearing clearance exaggerated). 
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The last matrix in Eg. (10) is known as the “stiffness matrix,” 
since it stands as a multiplier of the rotor displacement vector. 
The diagonal terms of this matrix are, indeed, the stiffnesses of 
the rotor and of corresponding fluid film in the x and y direc- 
tions. The off-diagonal terms in this matrix are known as cross- 
coupled stiffnesses,” and are responsible for rotor instability. As 
can be seen, these skew-symmetric terms are products of fluid 
damping D A, and rotative speed 9, each of them having a clear 
physical sense. The “cross-coupled stiffness” is simply a result 
of fluid viscous damping rotation at the circumferential velocity 
AQ, 


Fluid-force model with fluid inertia effect 

Consider now that the fluid force (2) contains, not only the 
spring and dashpot, but also an inertia element (Fig. 4): 

F.=M,@+DE+K,5 (11) 
where M, is the fluid inertia effect, and F, is the fluid force 
rotating at the rate Ag. 

Following the same reasoning used previously, the rotor/fluid 
system model in the stationary coordinates is as follows: 

Mz+D, 2+ Kz+M,(z-2jAnz-A7a"z) + DEAR) +K, z=0 





(12) 

or in the matrix form: 

| M+M, ols D:+D 2M, aAgl|x 

+ | 

0 M+M AY -2M,A2 D.+DIP 

K+K MA O7 AQD| |x 
~ =0 (13) 

-AQD K+K,-M,A7971 | y 


As it can be seen, some new fluid-inertia related terms 
appeared in the rotor/fluid system model equations (12) and 
(13). The damping matrix in Eq. (13) is now complemented by 
the skew-symmetric, “gyrascopic’’-like terms due to the Cori- 


ig 


fe 


Figure 4 
Fluid force model with fluid inertia effect. 
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olis acceleration of the rotating fluid. The stiffness matrix diago- 
nal terms are decreased by fluid inertia centrifugal force, since 
M,A Q" appears with the minus sign. This implies that the 
resulting radial “stiffness” may become negative at some high 
rotative speed ©. And it really does, which was shown during 
modal perturbation testing [6]. 

The question remains whether the fluid damping and fluid 
inertia rotate at the same rate AQ. Most probably, in a general 
case, these rates are different [7], that is, Eq. (12) should be 
rewritten as follows: 

Mz +D Z +K +M, (2-2jA,Q72-A2072) + D(Z-jAQz) +K,Z=0 

(14) 


where AQ is the angular velocity at which fluid inertia force 
rotates. The results of several experiments indicate that A,is up 
to 10% larger than A, especially at higher rotative speeds [7]. 


Final remarks 


The model of a rotor with fluid interaction (14) is still the 
simplest, as it assumes only one isotropic lateral mode of the 
rotor, and the fluid force action at the rotor modal mass loca- 
tion. The model can be extended to include more lateral modes 
of the rotor, and to consider rotor and fluid force anisotropy. 
Finally, the fluid force, which was assumed here to be linear 
with displacement and velocity, should be expanded to include 
its nonlinearity, especially the nonlinearity versus displacement. 


This article shows how to start building the rotor/Muid 
models. It emphasizes the physical meaning of involved param- 
eters, and treats the rotor and fluid as one system. The model in 
the form of a differential equation is only the first step. It has to 
be solved in order to trace the actual behavior of the rotor, its 
lateral displacement under the action of the fluid force. The 
solution of Eq. (14) will be discussed in a following article. 
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